Quantum Critical Phenomena of ""He in Nanoporous Media 
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The superfluid transition in liquid *He filled in Gelsil glass observed in recent experiments is discussed in the 
framework of quantum critical phenomena. We show that quantum fluctuations of phase are indeed important at 
the experimentally studied temperature range owing to the small pore size of Gelsil, in contrast to *He filled in 
previously studied porous media such as Vycor glass. As a consequence of an effective particle-hole symmetry, 
the quantum critical phenomena of the system are described by the 4D XY universality class, except at very low 
temperatures. The simple scaling agrees with the experimental data remarkably well. 



Introduction- Superfluidity is one of the most impressive 
macroscopic quantum phenomena that can be experimentafly 
observed. Since liquid "*He is a very clean system and the 
U(l) symmetry of the quantum phase is exact, its superfluid 
transition is also an ideal case to study phase transitions. In 
fact, much of the most precise experimental estimates of criti- 
cal exponents is obtained from the superfluid transition of "^He 

Qi. 

In bulk liquid "^He, the second-order superfluid phase tran- 
sition occurs at finite temperature only. On the other hand, 
the properties of liquid "^He can be changed by confining "^He 
in porous media. For example, '^He in Vycor glass exhibits a 
quantum phase transition as a function of film thickness, when 
"^He forms a two-dimensional film on the interior surface of 
the pores |3]. However, the physics might be complicated due 
to the crossover from two-dimensional to three-dimensional 
behaviors. When "^He is filled in Vycor pores, the superfluid 
transition remains at finite temperature (see [4] and references 
therein). Recently, however, by using Gelsil glass which has 
nanopores of 2.5 nm mean nominal diameter, considerably 
smaller than comparable substrates used in the past, the su- 
perfluid transition temperature is suppressed to zero under an 
applied pressure |5]. This implies the existence of a quantum 
critical point at zero temperature. 

The non-superfluid phase close to the quantum critical point 
exhibits rather abnormal properties. This also seems to be 
closely related to the pseudogap regime in high-Tc supercon- 
ductivity iH. The ^He in nanoporous media is in several ways 
(including the absence of fermionic degrees of freedom) sim- 
pler, and might provide a useful insight into the physics of 
more complex systems such as the high-Tc superconductors. 

In this Letter, we discuss the quantum critical phenomena 
of the superfluid transition in liquid "*He in Gelsil glass. We 
determine the value of the interaction strength V in the Bose- 
Hubbard model for liquid "^He in porous media. Although the 
strongly interacting and dense quantum liquid nature of '^He 
makes a microscopic theoretical calculation of the effective 
parameter V difficult, we show that it can be derived from 
macroscopic properties of liquid "^He. We find that, although 
the interaction parameter is about 20 mK for the Vycor sys- 
tem, it is as large as about 1 K for the Gelsil system. This 
explains the fact that the quantum critical phenomena are vis- 



ible for the experimentally studied temperature range only in 
the Gelsil system. 

The Bose-Hubbard model with disorder has been a subject 
of extensive theoretical study. Its quantum critical phenomena 
are highly nontrivial, and much of the problem still remains 
open. Nevertheless, we argue that the particle-hole symmetry 
breaking due to chemical potential, which drives the system 
away from the 4D XY universality class, is practically neg- 
ligible in the present system. As a consequence, for a large 
part of the phase diagram, the quantum critical phenomena 
observed experimentally in "*He in Gelsil glass can be under- 
stood with the simple non-random 4 dimensional XY (4D XY) 
universality class [7]. 

We demonstrate that the simple scaling based on the non- 
random 4D XY universality class indeed agrees quite well 
with the experimental data on liquid '^He in Gelsil glass, ex- 
cept in the small region close to the quantum critical point, 
where a crossover to a different universality class is observed. 

Setup- The Gelsil glass may be modeled by many highly 
interconnected pores, each of which can contain a number 
of "^He atoms. While the size of the pores is randomly dis- 
tributed, there is a typical pore diameter for a given sample. 

One of the most outstanding experimental findings on the 
system is a rounded peak in the specific heat at temperature 
Tb somewhat lower than the superfluid transition temperature 
T\ in the bulk. In Ref. |4] it was proposed that the specific 
heat peak should be interpreted as the formation of localized 
BEC (LBEC) within each pore. 

In fact, here we demonstrate that the peak temperature Tb 
can be indeed understood as a rounded lambda transition tem- 
perature in a finite-size system. Standard scaling theory gives 



AT = Tx-TBil) oc Cr^/'^^", 



(1) 



where I is the linear pore size, and 1/3/) ^ 0.67 is the coiTela- 
tion length exponent for the A transition. 

In two samples of nanoporous glasses with different pore 
sizes, Gelsil and Vycor, different values of AT were ob- 
served ijsiQl- Moreover, a similar shift was also observed 
in porous Gold samples with two different pore sizes fg]. The 
four values of AT are plotted as a function of the pore size I 
in Fig.[I] Here we use the effective pore size, obtained by sub- 
tracting the inert layer thickness from the nominal pore diam- 
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eter. Although the data were taken for different porous mate- 
rials at different pressures, the agreement with the scaling ([T]i 
is remarkable, v^j^, — 1/1.48 — 0.676 estimated from the fit 
is consistent with the known value v^j:, ~ 0.67. In addition to 
establishing the localized BEC picture, this result implies that 
liquid "^He in the nanopores does inherit bulk properties. 




■ Gelsin = 1.3nm 
^ Vycor / — 5 nm 
Gold ; = 23 iim 
▼ Gold I = 74 urn 
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FIG. 1 . (color online) The condensation temperature shift T\ — Tb 
as a function of the effective pore size, which supports the scaling Q. 

Since the number of atoms is finite, "^He in a single pore 
cannot be a superfluid in the true thermodynamic sense. In 
fact, the superfluid transition temperature is lower than Tb- 
However, below Tb, liquid "^He in each pore can have an 
approximately defined quantum phase which is a character- 
istic of a condensate. On the other hand, superfluidity as a 
macroscopic phenomenon only occurs when, by the hopping 
of atoms between pores, global phase coherence emerges over 
the entire system. This occurs at a temperature Tc which is 
lower than Tb ■ 

The system may be represented by the (disordered) Bose- 
Hubbard model 
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V, 



- ^ [tijala^ +h.c. ) , (2) 



where each site corresponds to a pore containing a local con- 
densate, connected to neighboring sites by random hopping 
tij . Every pore has a random chemical potential pi and a ran- 
dom charging energy Vi. The indices i,j label the pores and 
fii represents the number of "^He atoms in the i-th pore. The 
charging energy and the chemical potential can be combined 
as 



const.. 



(3) 



where = fM/Vi. The quantum phase di at site i can 
be regarded as canonically conjugate to the number operator 
Tii, namely they satisfy the canonical commutation relations 
[6j,nk] = iSjk, implying that the phase and atom number in 
each pore obey the uncertainty relation. The "charging en- 
ergy" Vi tends to fix Ui to fii, thus introducing quantum fluc- 
tuations of the phase 6i due to the uncertainty relations. 



Quantum Fluctuations- It is important to estimate the typi- 
cal value of the "charging energy" Vi. In the application of 
the Bose-Hubbard model to superconductors, Vi represents 
the charging energy due to Coulomb repulsion. In contrast, 
in the present case of neutral liquid ''He, there is no Coulomb 
repulsion. Nevertheless, putting an extra atom to, or remov- 
ing an atom from, a pore in the groundstate should lead to an 
increase in energy. This can be related to the finite compress- 
ibility of the liquid in the pore. In fact, we can estimate Vi as 
follows: 



1 

V, 



drii 



(4) 



where Vi is the (effective) volume of the pore, and v and k are 
the number density and the compressibility of liquid ''He, re- 
spectively. We approximate v and n of "'He inside the pores by 
their values in the bulk: v ^ 2.1x 10^^ m~^, n ^ 10"'' Pa"\ 

Then, assuming that a typical pore is a sphere with effective 
diameter 2i?eff =1-3 nm (subtracting the thickness 0.6 nm of 
the inert layer), we obtain an estimate for the typical value of 
Vi ^ 1.4K for ■'He in Gelsil Glass. In contrast, for Vycor 
Glass with effective diameter 5 nm, the same argument leads 
to Vi ~ 0.02K. 

If T ^ Vi, the quantum fluctuations of the phase may be 
neglected. In such a limit, although the phase 9i is a quantum- 
mechanical degree of freedom, it can be regarded as a classical 
variable. On the other hand, when T < Vi, quantum fluctu- 
ations of the phase become important. Thus, in the temper- 
ature range probed in the experiments, quantum fluctuations 
would be important for ''He filled in Gelsil glass but not in 
Vycor glass. This is consistent with the experimental results 
that quantum critical phenomena are observed only in Gelsil 
glass, when pores are filled by ''He. 

In the continuum limit, the (disordered) Bose-Hubbard 
model ^ is described by the Lagrangian density in 3 + 1 
dimensions: 
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- ^r[dr - 90 - 59ix)r^P + |(ro + Srix)M^ 
^]uo\M\ (5) 



where is a complex scalar field, g^, tq, and uq are constants, 
and Sg{x) and Sr{x) are position-dependent random variables 
distributed around zero 1 3 U 1 • 

When go = Sg = Sr = 0, this theory reduces to the stan- 
dard ?/;^ theory with the dynamical critical exponent z — 1, 
and the quantum critical point is described by the 4D XY uni- 
versality class. Since 4 is the upper critical dimension for the 
XY model, the critical exponents are given by the mean-field 
theory. For example, the correlation length exponent is given 
by i^iD = 1/2 [12]. This limit corresponds to, in terms of 
the original model (|2]l, the special case where there is no dis- 
order and Ui — fii/Vi's are exactly an integer. In particular, 
go = Sg = implies an exact particle-hole symmetry. 

In reality, fii is randomly distributed, breaking the particle- 
hole symmetry. Its effect can be classified into the overall 
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symmetry breaking go and the local symmetry breaking 5g, in 
the continuum theory (|5]l. go is generically non-vanishing and 
is a relevant perturbation to the 4D XY fixed point. In fact, 
even in a system without randomness, drives the system 
to a different critical behavior except at special multicritical 
points (tips of the Mott lobes). Nevertheless, here we argue 
that the effect of go is practically negligible for the present 
system, partly owing to the randomness. 

Let us assume that the distribution of the pore radius has 
standard deviation of, say, Ai? ~ 0.02nm, which would be 
rather an underestimation. This translates to the width 



47rj/i?eft- Ai? ~ 2.2 



(6) 



for the distribution of the n;. We assume that the follow a 
Gaussian distribution with average n^v and standard deviation 
An. Since the effect of the particle-hole symmetry breaking 
is a periodic function of n.i, it may be estimated by sin27r7l,i. 
The effective overall particle-hole symmetry breaking go then 
reads 



sin {2tix)- 



dx ~ sin (27rnav)e 



(7) 



/27rAn 

The first factor sin (27rnav) just represents the particle-hole 
symmetry breaking for the average value n av, which is gener- 
ically non-vanishing. The second factor e"^'^ shows the 
suppression of the symmetry breaking by the random distribu- 
tion. 

For the width given in eq. (|6l), the suppression is in fact 
about 10^'*'^. Thus, in the realistic temperature range, the 
overall particle-hole symmetry breaking go is negligible. That 
is, in the present system, each pore contains enough particles 
so that the distinction between the particle and hole becomes 
unimportant. A similar discussion was given for the asymp- 
totic low-energy behavior of disordered bosons |11]. Here 
the situation is somewhat different in that the smallness of go 
is not a consequence of the RG transformation, but is rather of 
microscopic origin. 

The random particle-hole symmetry breaking 6g is an irrel- 
evant perturbation to the 4D XY (Gaussian) fixed point ll ill . 
Although the presence of 5g is believed to be important to de- 
termine the eventual fate of the RG flow, it can be ignored in 
the neighborhood of the 4D XY fixed point. Thus tq and 5r 
are the only remaining relevant parameters around the 4D XY 
fixed point. This theory represents the critical behavior of a 
classical XY model in 4 dimensions. In the absence of the 
random 5r, the quantum phase transition belongs to the 4D 
XY universality class, ro corresponds to the temperature of 
the classical XY model, and is the control parameter for the 
quantum phase transition. It is thus identified with the pres- 
sure p of "^He in the system. 

The random term 5r corresponds to disorder in the classical 
XY model. We assume that the disorder is spatially uncorre- 
lated. However, in the mapping to the classical XY model, the 
disorder is completely correlated in the imaginary time direc- 
tion. This class of disorder turns out to be a relevant perturba- 
tion to the 4D XY fixed point, driving the system to a different 



fixed point f\3\ . This is in contrast to the effect of disorder 
on the finite temperature superfluid transition, where the un- 
corrected disorder is an irrelevant perturbation according to 
the Harris' criterion Hi, 151. 

The new fixed point induced by 5r, corresponding to the 
four-dimensional XY model with disorder correlated in the 
imaginary time direction, is sometimes called the random rod 
fixed point. At this random rod fixed point, the random part 
of the chemical potential 6g is believed to be relevant, even- 
tually driving the system to yet another (disordered boson) 
fixed point i 1 li] . However, if 5r is small, the system may be 
described by the 4D XY fixed point and the crossover away 
from it due to 5r, down to a certain temperature. 

We argue that this is indeed the case in the present system 
of "^He in Gelsil glass, concerning the experimentally studied 
temperature range. In fact, as we will demonstrate in the fol- 
lowing, a large part of the experimental data fits quite well 
with the simple 4D XY scaling. At very low temperatures, 
we see a crossover from the 4D XY behavior, which would be 
primarily due to the 5r perturbation. 
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FIG. 2. (color online) Frequency shift in the torsional oscillator 
experiment, which is proportional to the superfluid density, extrap- 
olated to zero temperature and as a function of the pressure. The 
experimental data, shown in squares, are taken from Ref 1 16]. The 
solid line is the 4D XY linear scaling {S). The agreement between 
the data and the scaling is remarkable except in the close vicinity of 
the quantum critical point. 

Superfluid density at zero temperature- Let us discuss a few 
consequences of the 4D XY critical behavior, which we com- 
pare with the experiments. The superfluid density scales with 
the exponent C, = (d + z — 2)v at zero temperature. For the 
4D XY universality class, d = 3, z = 1 and the correlation 
length exponent is v^^o = 1/2. Thus we find 



p.(p,T = 0)a(pc(0)-p)^^° =pM-P 



(8) 



where Pc(0) is the critical pressure at zero temperature. We 
show the experimental data for the superfluid density extrap- 
olated to zero temperature in Fig. IH and compare with the 
prediction (O. The excellent agreement of the data with the 
linear scaling strongly supports our proposal based on the 4D 
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XY universality class. For example, if we had z = 2 instead 
of 2; = 1, it would follow ps oc (pc(0) —pf'^'^ which is clearly 
inconsistent with the data. 

We note that the extrapolation of the linear scaling to zero 
temperature gives Pc(0) ^ 3.0 MPa instead of the 3.4 MPa 
which is the critical pressure directly obtained from the exper- 
imental data. We suppose that the non-random 4D XY scaling 
is valid for the effective critical pressure Pc^(O) ~ 3.0 MPa; 
this is different from the "true" critical pressure Pc(0), which 
is likely to be affected by the randomness. 
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FIG. 3. (color online) Phase diagram as obtained in fl^ and the fit 
to the experimental data. 

Finite temperatures- At finite temperature T, the length of 
the system in imaginary time direction takes the finite value 
1/T and we impose periodic boundary conditions in this di- 
rection, see for example fl?'] . The finite temperature problem 
is thus equivalent to the classical XY model on a hyperstrip 
geometry where the base is an infinite three dimensional sim- 
ple cubic lattice and the slab width is given by 1/T. The sys- 
tem at finite t emp erature is thus amenable to a finite size scal- 
ing analysis 11811 in the finite imaginary time direction. This 
allows us to make several predictions on the temperature de- 
pendence of the system. 

For example, the critical pressure Pc{T) at temperature T 
follows 

pf{0) - PciT) cx T^/"'^" = T^. (9) 

In Fig. Owe fitted a power-law to the experimental data of the 
phase boundary between the su^erfluid and the non-superfluid 
(LBEC) phases presented in 116]. We find an exponent of 2.13 
which is in good agreement with equation (|9]), while p"^^ (0) ^ 
3.2 MPa is slightly larger than that obtained in Fig.|2] 

Discussion- We have demonstrated that the experimental 
data on ^He in Gelsil glass agrees with 4D XY scaling rather 
well. This corresponds to the theory ^ with a small pertur- 
bation Sr, displaying clean 4D XY quantum critical behavior 
sufficiently far away from the quantum critical point. How- 



ever, close to the quantum critical point, the effect of Sr is en- 
hanced and crossover to a different universality class should 
occur In fact, by inspection of Fig. [3] the phase boundary at 
very low temperatures starts to deviate in a pronounced way 
from the 4D XY scaling ©. We also find that, in Fig.|2] the 
superfluid density at zero temperature also deviates from the 
4D XY scaling ([8]), close to the critical pressure Pc{Q)- These 
effects of the disorder are the subject of future investigations. 

We also note that, close to the quantum critical point, we 
find that the extraction of physical quantities (such as the su- 
perfluid density and the critical temperature) from the raw fre- 
quency data also becomes rather subtle. A more precise analy- 
sis would require a more sophisticated analysis of the raw data 
based on a better theoretical understanding. In any case, we 
believe that the analyses presented in this Letter demonstrate 
the basic validity of the proposed picture. 
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